Microscopic analysis of shot-noise suppression in nondegenerate diffusive conductors 
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We present a theoretical investigation of shot-noise suppression due to long-range Coulomb inter- 
action in nondegenerate diffusive conductors. Calculations make use of an ensemble Monte Carlo 
simulator self-consistently coupled with a one-dimensional Poisson solver. We analyze the noise in 
a lightly doped active region surrounded by two contacts acting as thermal reservoirs. By taking 
the doping of the injecting contacts and the applied voltage as variable parameters, the influence 
of elastic and inelastic scattering in the active region is investigated. The transition from ballistic 
to diffusive transport regimes under different contact injecting statistics is analyzed and discussed. 
Provided significant space-charge effects take place inside the active region, long-range Coulomb 
interaction is found to play an essential role in suppressing the shot noise at qlJ ^ ksT. In the 
elastic diffusive regime, momentum space dimensionality is found to modify the suppression factor 
7, which within numerical uncertainty takes values respectively of about 1/3, 1/2 and 0.7 in the 3D, 
2D and ID cases. In the inelastic diffusive regime, shot noise is suppressed to the thermal value. 
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I. INTRODUCTION 

Shot noise is caused by the randomness in the flux of 
carriers crossing the active region of a given device, anc 
associated with the discreteness of the electric charge.l 
At low frequency (small compared to the inverse of the 
transit time through the active region, / ^ I/'T't, but 
sufflciently high to avoid 1// contributions) the power 
spectral density of shot noise is given by 5'/(0) — ^2ql, 
where / is the dc current, q the electron charge, and 7 
the suppression factor. Uncorrelated carriers exhibiting 
Poissonian statistics are known to be characterized by 
a full shot-noise power (7 = 1). However, correlations 
between carriers can reduce the noise, leading to sup- 
pressed shot noise with 7 < 1. Several interactions and 
mechanisms can introduce correlations among carriers, 
thus giving rise to different levels of suppression p which 
can provide valuable information concerning the carrier 
kinetics inside the devices not available, from dc charac- 
teristics or low-frequency conductanccd 

In ballistic systems, like vacuum tubes, shot noise is 
known since the seminal work of Schottkyjj and well 
understood in terms of the Poissonian statistics of in- 
jected carriers. Within this model, shot noise has been 
investigated also in several nonuniform devices like Schot- 
tky diodes, p-n junctions, tunnel diodes, etc.EI In con- 
trast with these ballistic or quasi-ballistic structures, in 
macroscopic devices, where scattering mechanisms with 
phonons, impurities and other carriers determine the 
transport properties, shot noise is not usually detected 
and noise levels close to the thermal value are typi- 
cally measured (in the frequency range beyond 1// and 



generation-recombination contributions) . 

With the advent in recent years of mesoscopic conduc- 
tors, shot noise is receiving renewed attention. In partic- 
ular, the phenomenon of suppression, being a signature 
of correlations among particles, has emerged as a subject 
of relevant interest. The suppression has been predicted 
theoretically as a consequence of Pauli exclusion prin- 
ciple under strongly degenerate conditions in very dif- 
ferent situations. |tB| ballistic regime, shot noise is com- 
pletely suppresscdQla due to the non-fluctuating occupa- 
tion number of the incoming states. In a point contact, 
a peak in the iipise is predicted in between the conduc- 
tance plateaus.|j In symmetric double-barrier junctions a 
1/2 suppression factor Jias been theoretically explained 
by different authors. Q~tj In the case of elastic diffusive 
conductors, a 1/3 reduction of th£,ncjise has been calcu- 
lated for noninteracting electrons Jiil'E3L3 while in the case 
of strong electron-electmrpcattering the value of the sup- 
pression factor is \/3/4Jl3'E£l Finally, when the devices be- 
come macroscopic and inelastic processes are present, like 
scattering with phonnns. the noise is expected to reduce 
to the thermal value.EHa Remarkably, many ofj-tiiese pre- 
dictions have been experimentally conflrmed,ll2l"E2l thus 
opening new and interesting perspectives. Within this 
scenario, the understanding of the physical mechanisms 
originating shot noise and its suppression in mesoscopic 
conductors, and more generally in small-dimensional de- 
vices, is a field of growing interest. 

Most of the theoretical works carried out so far have 
considered degenerate conductors, where the Pauli ex- 
clusion principle plays a major role, and neglected long- 
range Coulomb interaction among carriers. The influ- 
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ence of this interaction is known to be relevant pto the 
noise reduction since the times of vacuum tubesj^j and 
its possible role in the case of mesosc(mic--.samplesE3 
has been repeatedly claimed by Landauer.EjEj However, 
only recently some works, always dealing with degener- 
ate conductors, include explicitly long-range Coulomb 
interactioncH and analyze its mfljuence on the high- 
frequency spectrum of shot noise.Ej ~lj A systematic anal- 
ysis of shot-noise suppression in nondegenerate conduc- 
tors with the inclusion of Coulomb correlations is still 
lacking. 

Coulomb interaction can affect noise in two main ways. 
On the one hand, being responsible for total current con- 
servation, it leads to local voltage fluctuations which try 
to preserve charge-peutrality and as such may influence 
the current noise.Firn On the other hand, being a repul- 
sive force, it tends to space electrons more regularly than 
a Poissonian statistics, thus, reducing the possible noise 
present in the current flux.ES This regulation of the elec- 
tron motion is particularly evident in the case of ballistic 
transport, which was investigated under nondegenerate 
conditions in previous works JlJcJ However, to our opin- 
ion, the role of Coulomb interaction in suppressing shot 
noise in the presence of scattering is still not well as- 
sessed. Here two issues are of main concern: (i) the de- 
termination of the reduction value under elastic diffusive 
conditions, and (ii) the understanding of the progressive 
disappearance of shot noise when passing from the meso- 
scopic to the macroscopic realm of conduction under the 
influence of inelastic processes. 

The aim of this paper is to shed new light on the 
above issues by investigating microscopically shot noise 
and its suppression in nondegenerate conductors in the 
presence of elastic and inelastic scattering and long-range 
Coulomb interaction. Our approach differs from those 
typically used to analyze noise in mesoscopic systems. 
Calculations are based on an ensemble Monte Carlo (MC) 
simulation self-consistently coupled with a Poisson solver 
(PS). Here the scattering mechanisms and the fluctua- 
tions of the self-consistent potential are intrinsically ac- 
counted for. In addition, the approach can analyze dif- 
ferent voltage-bias conditions ranging from thermal equi- 
librium to high electric fields necessary for shot noise 
to mamfest without the difficulties that other methods 
meet.Ell With this approach we investigate shot-noise 
suppression under the following conditions: (i) crossover 
between ballistic and diffusive transport regimes and dif- 
ferent carrier injecting statistics, (ii) diffusive transport 
regime under elastic or inelastic scattering, and (iii) 3D, 
2D and ID momentum space. A MC simulation of shot- 
noise suppression in a ID mesoscopic conductor has been 
performed in Ref. It considered a degenerate sam- 
ple in the presence of elastic and inelastic scattering 
where Pauli exclusion principle was taken into account 
but Coulomb interaction was neglected. 

Despite the fact that the physical system we analyze is 
nondegenerate, some of the results we get coincide with 
those obtained in degenerate systems by using other ap- 
proaches. In particular, it is specially surprising that 
we find the same universal 1/3 value for the shotjuoise 
suppression factor under elastic diffusive regime. cZI This 



1/3 value has been obtained by very different theoret=. 
ical approaches, going from quantum-phas e-Gph | er |entE3 
to semiclassical-diffusive degenerate models.tilEJcj The 
universa Uty o f this factor has been demo nstr ated both in 
quantumEacj and semiclassical contexts. EiJEil In all these 
cases degenerate conditions are assumed, and the noise 
reduction comes from the regulation of electron motion 
by the exclusion principle. However, it is not clear if the 
different approaches are equivalent, and the reappearance 
of the 1/3 factor could be just a njumerical coincidence, 
as critically asserted by LandauerE3 Here we show an- 
other unrelated context where the 1/3 suppression fac- 
tor appears, in which neither phase coherence nor Fermi 
statistics are present. In our case the origin of the effect 
is completely classical, and the correlation betwjeen elec- 
trons comes just from their Coulomb repulsion.cJ More- 
over, we show that the 1/3 value only appears when a 3D 
momentum space is considered. In the 2D or ID cases, 
different factors are obtained, though the physical mech- 
anism of suppression remains the same. To this purpose, 
an analytical theory which explains this dependence of 
the suppression factor on the dimension of momc*3jtum 
space has been recently developed by Beenakker.E3 Fi- 
nally, we will also illustrate the essential role played by 
Coulomb interaction for the suppression of shot noise. by 
inelastic scattering, as already stressed by BiittikerH3 

The outline of the paper is as follows. In Sec. II we 
discuss the physical model used for the structures un- 
der analysis, with special emphasis on the modeling of 
the contacts. We also provide the details of the MC 
simulation and noise calculations. Section III presents 
the results of shot-noise suppression with reference to: 
(i) the crossover between ballistic and diffusive transport 
regimes, (ii) the elastic and inelastic diffusive regime and, 
(iii) the role of momentum-space dimensionality. In Sec. 
IV the main conclusions and future trends are outlined. 



II. PHYSICAL MODEL 

In this section we present the details of the structures 
here analyzed, the models used in MC simulations and 
the procedures for noise calculation. 

A. Analyzed structure 

For our analysis we consider the simple structure 
shown in Fig. |[ It consists of a lightly doped active 
region of a semiconductor sample sandwiched between 
two heavily doped contacts (of the same semiconduc- 
tor) which act as thermal reservoirs and inject carriers 
into the active region. The sample is assumed to have 
a transversal size sufficiently thick to allow a ID elec- 
trostatic treatment in the x direction and neglect the 
effects of boundaries in y and z directions. The dop- 
ing of the contacts ric is taken to be much higher than 
that of the active region Nu. The carrier density at the 
contacts corresponds to their doping concentration; all 
impurities are assumed to be ionized at the temperature 
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T = 300 K considered here. The contacts are assumed 
to have no voltage drop inside and to remain always at 
thermal equilibrium. Accordingly, when a voltage U is 
applied to the structure, all the potential drop takes place 
inside the active region, between the positions x — Q and 
X — L. We shall analyze shot- noise suppression in differ- 
ent structures of same length but variable contact doping. 

B. Contact models 

The modeling of carrier injection from contacts can 
be crucial for the noise behavior in mesoscepic devices, 
especially in the case of ballistic transport .CJ To provide 
a complete model for the contacts and thus define the 
related sources of randomness in the carrier flux, we have 
to speciiy the velocity distribution of the injected carriers 
finji'v), the injection rate F and its statistical properties. 
Here we have denoted v = {vx,Vy,Vz)- 

Let us consider the process of electron injection from 
contact 1 into the active region at x—0 (see Fig. Q). Ac- 
cording to the equilibrium conditions of the contacts, 
the injected carriers follow a Maxwellian distribution 
weighted by the velocity component normal to the 
surface of the contact 

/mj(v) = Va:/Afs(v), > 0, (1) 

where /m b (v) is the Maxwell-Boltzmann distribution at 
the lattice temperature. The injection rate F, i.e. the 
number of carriers per unit time which enter the sample, 
is given by 

F = UcV+S, (2) 
where S is the cross-sectional area of the device, and 

.oo .oo .oo [k^ 

v+ = / / Jinj[^)avxdvydvz = \ , (3) 

with ks the Boltzmann constant and m the carrier effec- 
tive mass. The injection rate is taken to be independent 
of the applied voltage. Due to the very high value of ric 
as compared to iVu, any possible influence of the applied 
voltage (especially for high values) on the contacts and, 
consequently, on the injection rate, is neglected. Thus, 
the boundary injecting condition at the contacts is de- 
scribed through the constant injection rate F. The max- 
imum current that a contact can provide in the ballistic 
limit is the saturation current Is — qT. 

According to the nondegenerate distribution of carri- 
ers, the random injection at the contacts is taken to fol- 
low a Poissonian statistics. Thus, the time between the 
injection of two consecutive electrons, Unj, is generated 
with a probability per unit time given by 

P{Unj)=Te-^''"K (4) 

In the simulation we make use of Eq. (^) to generate 
tinj, which, following the MC technique, is given by 
tinj = where r is a random number uniformly 

distributed between and 1. Electrons are injected at 



a; = and x = L into the active region of the struc- 
ture according to the above stochastic rate. When a car- 
rier exits through any of the contacts it is canceled from 
the simulation statistics, which accounts only for carri- 
ers that are inside the active region at the given time 
t. Thus, the number of carriers in the sample N{t) is a 
stochastic quantity which fluctuates in time due to the 
random injection from the contacts and we can evaluate 
both the time-averaged value (N) and its fluctuations 
SN{t) = N{t)- {N). 

Consistently with the condition assumed here that the 
doping of the contacts is much higher than that of the 
active region, all the built-in effects associated with the 
diffusion of carriers around the contacts take place exclu- 
sively in the active region. Therefore, the effects related 
to possible charge fluctuations at the contacts are ne- 
glected in the calculation of the current. In any case, 
these effects are expected to appear at high frequen- 
cies (comparable with those of the plasma) while we are 
mostly interested in the low-frequency region of the noise 
spectrum. 

Unless otherwise indicated, calculations shall make use 
of the above contact model, which appears to be phys- 
ically plausible under nondegenerate conditions. How- 
ever, to analyze the influence of the contact injecting 
statistics on the noise behavior, alternative models are 
also used. In particular, for the injected carriers we shall 
consider: (i) fixed velocity instead of Maxwellian distri- 
bution and, (ii) uniform-in-time instead of Poissonian in- 
jection. In case (i) we consider the same injection rate F 
as in the basic model, but all carrier s are injected with 
identical cc- velocity Vx = a/ TrksT /2m^ which corresponds 
to the average velocity of the injected electrons when they 
follow a Maxwellian distribution. In case (ii) carriers are 
injected into the active region equally spaced in time at 
intervals of 1/F. 

C. Monte Carlo simulation 

The transport analysis is carried out by simulating the 
carrier dynamics only in the active region of the struc- 
ture. The influence of the contacts is included by means 
of the stochastic injection rate taking place at positions 
X — Q and x — L as described previously. Under the 
action of a dc applied voltage U , the carrier dynamics is 
simulated by an ensemble MC technique self-consistently 
coupled with a PSEJ The simulation is ID in real space, 
the Poisson equation being solved in the direction x of 
the applied voltage. Typically, a 3D momentum space is 
considered. However, to analyze the influence of dimen- 
sionality on the noise suppression, a 2D and ID momen- 
tum space is also considered in some specific cases. Since 
we are interested in analyzing only the effects related to 
Coulomb interaction, the electron gas inside the sample 
is assumed to be nondegenerate thus excluding any addi- 
tional correlations due to Fermi statistics. Carriers move 
inside the active region according to the semiclassical 
equations of motion with a constant effective mass, and 
undergo isotropic scattering in momentum space. Under 
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the condition of a constant appHed voltage, the instan- 
taMous current in a one-dimensional structure is given 
byH 

N(t) 

nt) = lY. «-(^)' (5) 

where Vxi{t) is the velocity component along the field 
direction of the «-th particle. It must be stressed that, 
although not explicitly appearing in Eq. (^) , the displace- 
ment current, which is of crucial importance for noise 
calculations, is imnlicitly taken into account by constant 
voltage conditions O 

To better analyze the importance of Coulomb correla- 
tions, which act by means of the fluctuations of the self- 
consistent potential, we provide the results for two differ- 
ent simulation schemes. The first one uses a dynamic PS, 
which means that the potential is self-consistently up- 
dated by solving the Poisson equation at each time step 
during the simulation to account for the fluctuations as- 
sociated with the long-range Coulomb interaction. The 
second scheme uses a static PS, thus calculating only the 
stationary potential profile; i.e., once the steady state is 
reached, the PS is switched off, so that carriers move in 
the frozen (non-fluctuating) electric field profile. Both 
schemes are checked to give exactly the same steady- 
state spatial distributions and total current, but the noise 
characteristics are different. Of course, the PS scheme 
which is physically correct is the dynamic one. The static 
scheme is just used to evaluate quantitatively the influ- 
ence of the potential fluctuations on the total noise, thus 
enabling us to separate the contribution belonging to ve- 
locity and number from that due to the self-consistent 
field. 

For the calculations we have used the following set of 
parameters: T — 300 K, m = 0.25rnoi dielectric constant 
e — 11.7£o (eo being the vacuum permittivity), L = 2000 
A, ric ranging between 10^'^ cm^'^ and 10^* cm~^ and 
Nd = 10^^ cm~^. According to Eq. (||), F is propor- 
tional to ric and it determines the level of space charge 
inside the active region, which will be c harac terized by 
the dimensionless parameter A, defined agl3c3 

A =7^, (6) 

^Dc 

where L^c = \/ eksT / q^ric is the Debye length corre- 
sponding to the carrier concentration at the contacts. In 
present calculations A will take the minimum value of 0.15 
{ric = 10^"^ cm^"^), for which the effects of Coulomb re- 
pulsion between electrons are practically negligible, and 
the maximum values of 30.9 {ric = 4 x 10^^ cm~^) and 
48.8 {ric — 10^^ cm^'^), for which quite significant elec- 
trostatic screening takes place. 

Scattering mechanisms are introduced in the simula- 
tion in a simple way by making use of an energy inde- 
pendent relaxation time r. Accordingly, we will consider 
separately elastic and inelastic interactions, both taken 
to be isotropic. In the elastic case, after each scattering 
event the direction of the velocity is randomized while 
the energy is conserved. In the inelastic case, beside 



randomizing its velocity, the electron is also completely 
thermalized by generating its velocity components after 
the scattering in accordance with a Maxwell-Boltzmann 
distribution at the lattice temperature. This model of 
strong inelastic scattering is adopted to evidence clearly 
the influence of energy dissipation on the noise suppres- 
sion. 

While L remains constant, the value of r is appropri- 
ately varied from 10 ps to 1 fs to cover both the ballistic 
and diffusive transport regimes. The transition between 
these regimes will be characterized by the ratio between 
the carrier mean free path ^, taken as H = '^v+t, and the 
sample length L. The conditions £/L ^ 1 and £/L <C 1 
correspond to the perfect ballistic and perfect diffusive 
regimes, respectively. 

Typical values of the time step and number of meshes 
in real space used for the PS are 2 fs and 100, respectively, 
except for the cases when r < 5 fs for which the time 
step is taken of 0.2 fs. As test of numerical reliability 
we have checked that by reducing the time step or by 
increasing the number of meshes the results remain the 
same. The average number of simulated particles in the 
active region ranges between 50 and 2000 depending on 
contact doping, transport regime, and applied voltage. 

D. Noise calculation 

In the analysis of shot noise we are particularly inter- 
ested in the low-frequency value of the spectral density 
of current fluctuations 5/(0), for then obtaining the sup- 
pression factor 7 = Si{0)/2ql. To this end, from the 
simulation we firstly evaluate the autocorrelation func- 
tion of current fluctuations 

Ci{t) - {SI{0)SI{t)), (7) 

where SI{t) — I{t) — (/) is the instantaneous current 
fluctuation. Then, the corresponding spectral density is 
obtained by Fourier transforming 

/oo 
Ci{t)e'^^f'dt. (8) 
-oo 

Once steady-state conditions in the active region of the 
structure are reached, the simulation typically consists 
of 350000 time steps. At the end of each time step the 
instantaneous current I{t) is calculated. From these suf- 
flciently long sequence of values of I{t), the time average 
of the current (/) is determined and the current autocor- 
relation function is directly calculated following Eq. (|^). 
To clarify the role of number and velocity contributions 
to the current noise, the autocorrelation function and the 
spectral density are decomposed into three main terms as 

Ci (t) = Cv {t) -f Cn (t) + CvN (t) , (9) 

respectively given by 

Cv{t) - ^{Nf{5v{Q)5v{t)), (10a) 

CN{t) = ^{vf{5N{Q)5N{t)), (10b) 
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CvN{t) = ^{v){N){SviO)SN{t) 

+ dN{0)5v{t)). (10c) 

In the above equations Cy is associated with fluctuations 
in the mean carrier velocity v(t) = J2^i Cn 
with fluctuations in the car ri a r pumber N(t), and Cvn 
with their cross correlationSEJ To distinguish between 
the results obtained from the static and dynamic PS, we 
shall denote the corresponding current spectral densities 
as Sf and Sf, respectively. 

III. RESULTS 

The main results of the present work are organized as 
follows. Subsections A and B pertain to a 3D momen- 
tum space while subsection C is devoted to 2D and ID 
momentum spaces. Most of the reported results refer to 
high values of the space-charge parameter A (typically 
A=30.9), which implies significant effects related to long- 
range Coulomtwjiteraction present in the active region of 
the structuresE3 The essential influence of A on the sup- 
pression factor 7 is analyzed at the end of subsection B. 
We recall that elastic and inelastic scattering are consid- 
ered separately in the simulations. In no case both types 
of scattering are taken into account simultaneously. 

A. Transition ballistic-diffusive regimes 

The behavior of noise in the crossover from ballistic to 
diffusive transport regimes is analyzed under far-from- 
equilibrium conditions {U ^ kBT/q), since these are 
necessary for the manifestation of shot noise. The depen- 
dence of noise on applied voltage for the perfect ballistic 
regime was already reported in Rcf. ^ and for the perfect 
diffusive regime it will be analyzed in next subsection. 

Fig. ^ reports the variance of carrier-number fluc- 
tuations normalized to the average carrier number 
(AA^2)/(Ar) as a function of i/L. The values taken by 
the variance are found to be practically independent of 
the transport regime. The results of the elastic and in- 
elastic cases are very similar; by contrast two different 
values are obtained depending on the scheme of the PS 
used in the calculations. When the static PS is consii. 
ered, it is found that, within numerical uncertainty,E3 
{AN'^) = (N) in aU the range of £/L reported. There- 
fore, in the absence of dynamic Coulomb correlations, 
the carrier number follows a Poissonian statistics. In the 
perfect ballistic regime this is a consequence of the in- 
jecting statistics at the contacts. In the diffusive regime 
this is due to the effect of the randomness introduced by 
scattering, independently of its elastic or inelastic prop- 
erty. Because of the Poissonian statistics observed for the 
carrier-number fluctuations, no shot-noise suppression is 
expected within the static PS scheme. On the contrary, 
when the dynamic PS scheme is used, a sub-Poissonian 
behavior is clearly observed with (AA^^) being about a 
factor of 0.4 lower than (A^); thus shot- noise suppression 
is expected. 



Fig. 1^ reports SijO) normalized to 2qls under the same 
conditions of Fig. ||. Here the evolution of the current in 
terms of 2ql is also shown. This evolution exhibits two 
limiting behaviors: saturation at £/L > 10~^, typical of 
a ballistic or quasi-ballistic regime; and linear decrease at 
£/L < 10^^ as £/L 0, typical of a diffusive regime. In 
both the elastic and inelastic cases, Si{0) calculated with 
the static PS coincides exactly with 2ql,ti thus leading 
to the conclusion that the stream of uncorrelated carriers, 
even in the presense of intensive scattering processes (no 
matter elastic or inelastic), exhibits full shot noise behav- 
ior like in the pure ballistic transport regime. On the con- 
trary, with the dynamic PS >5'/(0) is systematically lower 
than 2ql, thus evidencing a suppression effect. Here, as 
expected, elastic and inelastic curves recover the same 
value in the ballistic limit, the suppression correspond- 
ing to that caused by the fluctuations of the potential 
barrier near the cathode (induced by the ^p^fc, charge) 
which controls the current in this regime.cTE3 As the 
perfect diffusive regime is approached, the suppression 
remains active, more pronounced in the inelastic case, 
and is related to the joint action of Coulomb repulsion 
and the presence of scattering. 

We shall now consider the dynamic PS scheme only. 
Fig. ^ reports 7 as a function of £/L for several values 
of the applied voltage. In the perfect ballistic regime 
the two distinct values of 7 found refer to the presence oil 
absence of the potential barrier related to space charge.H 
For U = AOksT/q the barrier is still present and the 
suppression is important. For U = 80 and lOOfcsT/q the 
barrier has already disappeared; accordingly the current 
saturates and the suppression factor takes on the full 
shot-noise value. When the diffusive regime is achieved, 
in the elastic case 7 attains a constant value at further 
decreasing of £/L, and takes the same value of about 1/3 
independently of the applied voltage. On the contrary, 
in the inelastic case the higher the applied voltage the 
lower the value which 7 is found to take. Remarkably, 
the value of £/L at which 7 starts decreasing when the 
ballistic regime is abandoned is the same in the elastic 
and inelastic cases for a given applied voltage {£/L « 0.3 
and 0.1, for 80 and 100 ksT/q, respectively). However, 
when the diffusive regime is approached, a lower value 
of £/L is required in the inelastic case with respect to 
the elastic case for 7 to reach a constant value. This 
behavior can be explained in terms of the different elastic 
and inelastic scattering intensity required by the electron 
system to achieve a significant equipartition of energy 
into the three directions of momentum space. 

To better illustrate the above features. Fig. || shows 
the spatial profiles of the total average energy and its 
three contributions in the x, y and z directions of mo- 
mentum space for several values of t in the elastic and 
inelastic cases at C/ = lOQkBT/q . For the longest r, the 
total energy practically coincides with the energy in the 
X direction independently of elastic and inelastic cases, 
since the transport is nearly ballistic. In the case of elas- 
tic scattering [Fig. ||(a)], as t is shortened the energy in 
the X direction decreases while that in the other two di- 
rections increases, in such a way that the total energy 
remains constant, as follows from elastic conditions. By 
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further shortening r, the energy in the three directions 
becomes finally the same and equipartition condition is 
achieved. In the inelastic case [Fig. ^(b)], equipartition 
of energy is also reached, but only at the shortest r. The 
fact that energy equipartition is more easily reached in 
the elastic than in the inelastic case is understood as fol- 
lows. The occurrence of an elastic scattering provides to 
the y and z directions a part of the energy that the car- 
rier initially had in the x direction, so that the energy 
tends to equalize in the three momentum components. 
However, for the case of an inelastic scattering only the 
energy in the x direction changes, since the energy in the 
y and z directions always correspond to that of the ther- 
mal equilibrium distribution {^UbT). This means that 
for a given applied voltage, a higher number of scattering 
events (i.e. a lower value of l/L) is necessary to achieve 
energy equipartition in the inelastic case with respect to 
the elastic one. Thus, in Fig. |(a) it can be observed 
that for i/L = 5.38 x 10"'^ and an applied voltage of 
lOO/csT/g, in the elastic case an important isotropic re- 
distribution of energy has already been achieved, so that 
7 already takes the value of 1/3 [Fig. ^. However, in 
the inelastic case [Fig. 1(b)], for the same value of i/L 
the energy in the x direction is still significantly higher 
than that in the other two directions of momentum space 
and, as a result, 7 has not still reached its constant value 
(which corresponds to thermal noise, as will be shown 
later) and continues decreasing. These results confirm 
that thermal conditions, and therefore thermal equilib- 
rium noise, are not reached until the inelastic scattering 
time is so short that the energy gained by electrons in a 
mean free path is much lower thaii. the thermal energy, 
as already indicated by LandauerEj 

To analyze the role played by the modeling of the con- 
tact injection on the suppression of noise. Fig. ^ reports 7 
as a function of l/L calculated using four different mod- 
els. They combine Poissonian/uniform injecting statis- 
tics and Maxwellian/fixed- velocity distribution of the in- 
jected carriers. The Poissonian-Maxwellian injection is 
the basic one used in calculations. In the perfect bal- 
listic regime, when carrier transport is deterministic, 7 
crucially depends on the injection model. Thus, in the 
case of the uniform fixed-velocity model, when the injec- 
tion introduces no extra noise in the current flux, 7 is 
found to decrease linearly with the increase of £/L. The 
noise does not vanish completely since, unless l/L 00, 
there is always some probability of undergoing a scatter- 
ing mechanism. In this limit, when the noise is produced 
just by a few scattering events, it is clearly observed that 
elastic interactions lead to more important current fluc- 
tuations than inelastic mechanisms. By approaching the 
perfect diffusive regime the suppression factor is found to 
be independent of the model used. We conclude that the 
noise in the diffusive regime (and ■particularly the 1/3 
suppression value obtained in the elastic diffusive case) 
is independent of the carrier injecting statistics^ and it is 
only determined by the joint action of scattering mecha- 
nisms and Coulomb correlations. 



B. Diffusive regime 

In this section we shall consider scattering times short 
enough to ensure a diffusive transport regime {^/L < 
3 X 10""^). In this regime the noise behavior is clospfc, 
related to the breadth of the velocity distributionjlaa 
as shown for the x-velocity component in Fig. ^ In the 
case of elastic scattering, the distribution broadens at 
increasing applied voltages since there is no energy dis- 
sipation. In the inelastic case, the energy is maximally 
dissipated by reducing it to the lattice value after each 
scattering event; thus a thermal equilibrium Maxwell- 
Boltzmann distribution is obtained independently of the 
applied voltage. In both cases the distributions are very 
slightly displaced to positive velocities, as implied by the 
presence of a net current flowing through the structure. 
We remark that the distributions shown in Fig. ^ refer to 
all the carriers present inside the sample. In the elastic 
case, the local distributions at given positions are found 
to exhibit a nearly isotropic but no longer Maxwellian 
shape, however their spatial integration over the whole 
sample length gives the Maxwellian profile reported in 
the figure. 

Fig. H shows 5/(0) normalized to 2qls as a function 
of the applied voltage. Here, the I — U characteristic is 
also shown in terms of Sgl— In the wide range of volt- 
ages reported in the figure)^ the I — U curve exhibits a 
super-linear behavior which is related to the importance 
of space-charge effects and the increase of the carrier 
number inside the active region with the applied volt- 
age for qll > ksT. For the highest applied voltages it is 
found I oc Lf , where the power r is a function of A. In 
particular, r — 1.7 for A — 30.9 (the case shown in the 
figure) and r = 1.8 for A = 48.8. In what concerns S'/(0), 
near thermal equilibrium conditions {qU < ksT), elastic 
and inelastic cases exhibit the same value which satisfies 
Nyquist relation. In the inelastic case, at increasing volt- 
ages 5/(0) remains practically constant, in accordance 
with the velocity distributions of Fig. 0(b). The slight 
increase shown at the highest voltages is strictly related 
to the increase of the carrier number. On the contrary, 
in the elastic case at high voltages {qU /ksT > 10), when 
the velocity distribution starts broadening significantly 
from its thermal equilibrium shape [Fig. [7|(a)], 5/(0) in- 
creases systematically with J7, its ratio to the current 
remaining constant and providing a value of 1/3 for the 
suppression factor 7. 

The suppression factor as a function of the applied 
voltage is shown in Fig. ^(a) for A = 30.9, which cor- 
responds to the same conditions of the previous figure, 
and in Fig. ^(b) for A = 48.8. The behavior of 7 is quite 
similar for both values of A. By comparing the results 
corresponding to elastic and inelastic scattering we find 
that at the lowest bias both cases coincide, by providing 
the standard thermal noise given by Nyquist relation. 
On the contrary, at the highest voltages the elastic case 
reaches the 1/3 limiting value, while the inelastic case 
decreases systematically. In this latter case the values of 
7 are closely fitted by this simple expression for 5}"'^'(0): 
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where 



Go 



(11) 



(12) 



is the conductance and (iV)o the average number of elec- 
trons—inside the sample, both in the limit of vanishing 
biasO Thus, in the case of inelastic scattering the spec- 
tral density can be expressed analogously to that of ther- 
mal Nyquist noise modulated by the variation of carrier 
number at the given voltage, even in the presence of a 
high bias and a net current flowing through the struc- 
ture. We conclude that inelastic scattering strongly sup- 
presses shot noise and makes the noise become macro- 
scopic (7 <^ 1). These results, confirm previous prie-, 
dictions by Shjimizu and Ueda,ll3 Liu and Yamamoto,ta 
and Nagaev,Ej with the important difference that here 
we consider a nondegenerate classical conductor. We re- 
mark that present findings prove also that the condition 
of inelastic scattering alone does not suffice to suppress 
shot noise; the presence of the fluctuating self-consistent 
electric field remaining a necessary condition. Indeed, 
as a counter-proof we refer to the calculations performed 
with the static PS scheme, where no suppression has bew* 
detected (see Fig. ||). Therefore, as argued by Buttiker,EHl 
it is the combination of both Coulomb interaction and in- 
elastic scattering which leads to the suppression of shot 
noise. 

In the elastic case, the values of 7 are nicely reproduced 
by the following expression for Sf{Q): 



Sf{0) 



— K nTGn -y — r~ 

3 °{N)„ 



— ql coth I 

3^ \2kBT 



(13) 



which is quite similar to that obtained by Nagaevtu m 
a degenerate context. Eq. |l^ describes the crossover 
from thermal-Nyquist noise for qU ^ fc^T (where 
Sf = AkBTGo) to suppressed shot noise for qU ^ 
ksT (where Sf = §9-^)' l^us providing 7 — 1/3 
for the highest anpligd . v oltages. In contrast with 
other approaches ,E§ti 11 c3 our results show that, nei- 
ther phase- coherenceE^Ej nor degenerate statistics are re- 
quired for the occurrence of suppressed shot noise in dif- 
fusive conductors, and purely classical physical processes 
can lead to the same 1/3 factor 

To illustrate the physical origin of the 1/3 value, 
Fig. 0(a) reports a typical spectrum of the suppression 
factor under elastic diffusive conditions for static and dy- 
namic PS schemes. Here the current spectrum is decom- 
posed into velocity, number, and cross-correlation contri- 
butions, Si{f) = Svif) + SN{f) + SvNif), according to 
Eqs. ([lO|). In the static PS scheme the spectrum clearly 
shows that all the three terms contribute to 5/(/), and 
two different time scales can be identified. The longest 
one is associated with the transit time of carriers through 
the active region « 5 ps, and is evidenced in the terms 
Sffif) and 5'yjv(/). The shortest one is related to the 
relaxation time of elastic scattering t = 5 fs, and is man- 
ifested in Sy{f). Remarkably, the velocity contribution 



yields 1/3 of the full shot-noise value, while the others 
two terms provide the remaining 2/3. Thus, in the static 
PS scheme full shot noise is recovered as sum of all the 
three contributions. On the contrary, in the dynamic PS 
scheme 5^(/) and S'yjY(/) are found to exactly compen- 
sate each other and, as a result, Sf{f) coincides with 
S'y(/) in all the frequency range. Moreover, Sff{f) takes 
values much smaller than Sf^{f). The characteristic time 
scale of Sff{f) and Syj^{f) differs from that found within 
the static PS scheme, which was related to the transit 
time tt- Now in the dynamic case it is the dielectric re- 
laxation time corresponding to the carrier concentration 
at the contacts Td = 0.46 ps which determines the cutoff 
of the contributions belonging to number fluctuations. In 
the frequency range between the transit and collision fre- 
quency values it is interesting to notice that both static 
and dynamic PS schemes yield 7 — 1/3, thus relating 
the suppression factor to velocity fluctuations only. How- 
ever, at low frequencies only the dynamic scheme takes 
this value by virtue of Coulomb correlations, which are 
responsible for the reduction of S'^(/) and the mutual 
compensation of Sfi{f) and S'y^(/) contributions. It 
is remarkable that S'y(/) = Sy{f) in all the frequency 
range, which implies that velocity fluctuations are not 
affected by long-range Coulomb interaction, but just by 
scattering mechanisms. Coulomb repulsion affects only 
the contributions where carrier-number fluctuations are 
involved. Fig. 0(b) reports the spectrum for inelastic 
scattering. Here, the same features of the elastic case 
are observed, with the important difference that Sy{f) 
is much lower than when there is no energy dissipation. 

So far we have analyzed structures where Coulomb re- 
pulsion plays an important role (i.e. A ^ 1). To check 
to which extent this interaction is determinant for noise 
suppression. Fig. |l^ reports 7 and the three contribu- 
tions in which it has been decomposed as a function of 
A. Here, we present both the elastic and inelastic cases 
calculated under far-from-equilibrium conditions within 
the dynamic PS scheme. The variance of the number of 
carriers inside the sample is also shown to analyze the 
evolution of carrier statistics as a function of A. For the 
lowest values of A, when space-charge effects and in turn 
Coulomb interaction are negligible, full shot noise is ob- 
served. As A increases, 7 starts decreasing from unity 
until reaching a constant value for A 30. It is remark- 
able that the contribution of velocity fluctuations to 7 
does not vary signiflcantly with A, being much smaller 
in the inelastic case with respect to the elastic one be- 
cause of the quasi-thermal conditions imposed by energy 
dissipation. On the contrary, the contributions associ- 
ated with number and velocity-number fluctuations are 
strongly affected by the increase of A. Indeed, their ab- 
solute value decreases systematically and, being opposite 
in sign, they compensate each other at the highest values 
of A, so that Si = Sy. We also notice that the transition 
from full shot noise to suppressed shot noise regimes is 
sharper in the inelastic case. The variance of the carrier 
number inside the sample resembles the behavior of 7, in- 
dicating that the increase of A leads to more pronounced 
sub-Poissonian statistics until reaching the flnal regime 
for A ^ 1. Therefore, we conclude that the 1/3 value 
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exhibits several universal properties, namely, it is inde- 
pendent of: (i) the scattering strength (once £ L); (ii) 
the applied voltage (once qU :s> ksT), (ii) the screen- 
ing length (once A 1), and (iv) the carrier injecting 
statistics. 



C. Dependence on momentum space dimensionality 

The results reported so far refer to a 3D momentum 
space. In contrast to degenerate diffusive systems where, 
provided quasi-one dimensional conditions in real space 
are attained, noise suppression is independent of the 
number d of momentum space dimensions, an interesting 
feature of nondegenerate diffusive systems is that noise 
suppression can depend on d. For the inelastic case con- 
sidered here no dependence on d has been found, since 
there is no influence of the velocity components transver- 
sal to the electric field direction on transport and noise 
properties of the structures. On the contrary, in the elas- 
tic case the suppression factor is found to depend signif- 
icantly on df^ since the transversal velocity components 
constitute a channel for energy redistribution which af- 
fects the transport properties of the structure. Therefore, 
below we shall focus our analysis on the elastic case. Ac- 
cordingly, when d = 2 in the simulation the carrier veloc- 
ity is randomized in two components after each scattering 
event, and when d = 1 the isotropic character of scatter- 
ing is accomplished by inverting the carrier velocity with 
an average (back-scattering) probability Pb — 0.5. 

Fig. ^(a) reports 7 as a function of i/L for the ID, 2D 
and 3D cases at high voltages {U = iOksT/q) calculated 
within the dynamic PS scheme. We notice that, when 
calculated within the static PS scheme, the results in ID 
and 2D cases do not exhibit any shot-noise suppression, 
like in the 3D case. For the highest values oi £/L, in 
all three cases 7 approaches the asymptotic-value corre- 
sponding to the ballistic limit (7 = 0.045)c3, where the 
behavior is independent of d. At a given value of i/L, 
a higher deviation from the asymptotic ballistic value is 
observed for lower d. This is due to the fact that, in av- 
erage, elastic interactions introduce higher fluctuations 
of the carrier x- velocity the lower is the number of avail- 
able momentum states after the scattering mechanism (in 
particular, just 2 in the ID case). For this same reason, 
the increasing presence of scattering as i/L is reduced 
leads to higher values of the suppression factor the lower 
is the dimensionality. Remarkably, within numerical un- 
certainty the limit value reached by 7 in the perfect diffu- 
sive regime is found to be, respectively, of 1/3, 1/2, and 
0.7 for 3D, 2D, and 1D.E3 Fig. |l|(b), by reporting 7 in 
the diffusive regime as a function of the applied voltage, 
provides evidence that these limit values are independent 
of the bias once qU ^ ksT. The origin of the suppres- 
sion is the same in all three cases: the joint action of 
Coulomb correlations and elastic scattering, which leads 
to the result 5/(0) = Sv{0) [as shown in Fig. |l^(a) in the 
3D case], where Sv{0)/'2ql under perfect diffusive regime 
is a function of the dimensionality of momentum space. 

In a recent work, BeenakkerEj has developed an ana- 
lytical theory able to explain the dependence of 7 on d 



in nondegenerate diffusive conductors and has provided 
a closed expression for 7 as a function of d. In particular, 
the theory predicts for 7 the values of 0.34, 0.51, and 0.92 
for, respectively, the 3D, 2D and ID cases, in agreement 
(more closely for 3D and 2D) with present findings. From 
that analysis, Beenakker concluded that the proximity of 
7 to 1/d is accidental. In view of the simplifications in- 
troduced in his theory we believe that such a conclusion 
cannot rigorously be applied to the present model. 

For reason of completeness, we have finally analyzed 
the possible influence of an anisotropic elastic scattering 
on transport and noise in the ID case. Fig. |l^ reports 7 
and the current as a function oi i/L for three different 
values oi Pb'- 0.5 (isotropic case), 0.25 and 0.10 (less pro- 
nounced back-scattering). The noise results differ only 
in the quasi-ballistic regime, where the current remains 
nearly constant. Once the diffusive regime is achieved, 
the current decreases linearly with £/L and 7 takes the 
same value of about 0.7 independently of the degree of 
scattering anisotropy. As expected, we have found that 
the smaller the value of Pb the wider the quasi-ballistic 
range, and for a given value of i/L in this range the 
stronger the suppression, as corresponds to closer ballis- 
tic conditions. 



IV. CONCLUSIONS 

We have provided a microscopic analysis of shot-noise 
suppression in nondegenerate diffusive conductors. To 
this purpose, the carrier dynamics in the active region 
of a semiconductor structure under the influence of elas- 
tic or inelastic scattering has been simulated by using an 
ensemble MC technique self-consistently coupled with a 
PS. The essential role played by long-range Coulomb in- 
teraction on the shot-noise suppression factor 7 has been 
demonstrated, since no suppression is found in the ab- 
sence of the self-consistent potential fluctuations. 

We have analyzed shot-noise suppression in the region 
of crossover from ballistic to diffusive transport regimes. 
In the diffusive regime a value of 7 independent of sam- 
ple length is not achieved until a significant energy redis- 
tribution among momentum directions takes place. For 
high voltages qU/ksT ^ 1 and long samples £/L <C 1, 
in the elastic case shot noise is found to be suppressed to 
a 1/3 value, while in the inelastic case we have found a 
stronger suppression the higher the applied voltage. No- 
ticeably, in the perfect diffusive regime 7 is found to be 
independent of carrier injecting statistics, which implies 
that in this regime the noise is just a property of the 
sample. 

Our results show that neither phase coherence nor 
Fermi statistics are necessary for the appearance of the 
1/3 suppression factor in an elastic diffusive conductor. 
In our model, the appearance of this factor requires the 
simultaneous fulfillment of the three following conditions: 
e/L < 1, A > 1, and qU/ksT » 1. The first implies 
perfect diffusive regime, the second strong space-charge 
effects, and the third very far-from-equilibrium condi- 
tions. Inelastic scattering is found to further contribute 
in suppressing shot noise, by reducing it to values close 
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to thermal Nyquist noise under strong dissipative con- 
ditions. However, for this suppression to take place it 
is necessary the presence of long-range Coulomb interac- 
tion. 

The action of Coulomb repulsion in suppressing shot 
noise takes place through the reduction of the contribu- 
tions associated with carrier-number fluctuations to the 
total noise spectral density. In particular, the compensa- 
tion between number and velocity-number terms implies 
that the total noise is finally determined just by the con- 
tribution of velocity fluctuations. 

In the elastic case, 7 depends on the momentum space 
dimensionality, the suppression being less pronounced 
the lower the dimension of momentum space. 

We believe that the present investigation can provide 
a stimulus for an experimental verification of some of the 
reported results. As a valuable improvement we aim at 
including degenerate statistics in a further step. 
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FIG. 1. Schematic drawing of the structure under investi- 
gation. 



FIG. 4. Shot-noise suppression factor vs ballistic parameter 
£/L for the cases of elastic and inelastic scattering at different 
applied voltages. Calculations are performed by using the 
dynamic PS scheme. 



FIG. 5. Spatial profiles along the sample of the total aver- 
age energy and its three contributions in the x, y and z direc- 
tions of momentum space for several values of the scattering 
time r (and of the associated i/L). The applied voltage is 
U — lOOfcsT/g. (a) refers to elastic scattering, (b) to inelas- 
tic scattering. 



FIG. 6. Shot-noise suppression factor vs the ballistic pa- 
rameter £/L for an applied bias of U — AQknT/q calculated 
with different contact models. Calculations refer to the dy- 
namic PS scheme considering elastic and inelastic scattering. 



FIG. 7. Velocity distribution function of carriers inside 
the active region of the structure for l/L — 1.07 x 10~^ at 
several applied voltages, (a) refers to elastic scattering, (b) 
to inelastic scattering. 



FIG. 8. Low-frequency spectral density of current 
fiuctuations normalized to 2qls vs applied voltage for 
l/L = 1.07 X 10~^. Calculations refer to the dynamic PS 
scheme considering elastic and inelastic scattering. The cur- 
rent in terms of 2ql is also plotted for comparison. 



FIG. 9. Shot-noise suppression factor vs applied bias cal- 
culated with the dynamic PS scheme considering elastic and 
inelastic scattering with £/L — 1.07 x 10"'^ at two different 
levels of space charge: (a) A = 30.9 and (b) A = 48.8. The 
lines correspond to the fittings of Eqs. ( pT| ) and (p^). 



FIG. 10. Spectrum of the shot-noise suppression factor un- 
der diffusive regime (l/L = 2.69 X 10"^) calculated within 
static and dynamic PS schemes for (a) elastic and (b) inelas- 
tic cases and an applied voltage oi U = lOOkBT/q. Different 
contributions to the total spectrum are reported in the figure. 



FIG. 2. Variance of carrier number inside the active region 
normalized to the average carrier number vs the ballistic pa- 
rameter i/L for an applied voltage of !7 = AQknT /q. Calcula- 
tions are performed by using static and dynamic PS schemes, 
and considering elastic and inelastic scattering mechanisms. 



FIG. 11. Shot-noise suppression factor and variance of the 
carrier number inside the active region as a function of the 
characteristic parameter of space charge A under diffusive 
regime (£/L = 2.69 x 10""^) calculated within dynamic PS 
scheme for (a) elastic and (b) inelastic scattering. Velocity, 
number, and velocity- number contributions to the suppres- 
sion factor are also shown in the figure 



FIG. 3. Low-frequency spectral density of current fluctua- 
tions normalized to 2qls vs the ballistic parameter i/L for an 
applied voltage of C/ = AQksT /q. Calculations are performed 
by using static and dynamic PS schemes, and considering elas- 
tic and inelastic scattering mechanisms. The dependence of 
the current through 2ql is also plotted for comparison. 



FIG. 12. Shot-noise suppression factor for the cases of 1, 
2 and 3 dimensions of momentum space calculated within 
the dynamic PS scheme for elastic scattering as a function 
of: (a) ballistic parameter l/L with an applied voltage of 
U — lOfcsT/g and (b) applied bias U under diffusive regime 
{l/L = 1.07 X 10^^). 
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FIG. 13. Shot-noise suppression factor and normalized cur- 
rent vs the ballistic parameter tjL calculated within the dy- 
namic PS scheme for clastic scattering in the case of a ID 
momentum space. Different curves refer to the reported val- 
ues of the probability of scattering in the backward direction 
Pb- The appUed voltage is f/ = AQknT/q. 
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